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We study the nonequilibrium spin-thermal transport across a metal-magnetic insulator interface. The 
conjugate-converted thermal-spin transport is assisted by the exchange interaction at the interface, between con- 
duction electrons in the metal lead and localized spins in the insulating magnet lead. We predict the rectification 
and negative differential spin Seebeck effect and resolve their microscopic mechanism. The rectification of spin 
Peltier effect is also identified. The phenomena predicted here are relevant for designing efficient spin/magnon 
diode and transistor, which could play crucial roles in controlling energy and information in functional devices. 
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Recently, Spin Seebeck effect (SSE), a phenomenon that 
temperature bias can produce a spin current and an associ- 
ated spin voltage, has been observed in magnetic metals [1], 
magnetic semiconductors [2, 3], magnetic insulators [4, 5] and 
non-magnetic materials with spin-orbit coupling [6]. Since 
then, the SSE has ignited a new upsurge of research, because it 
acts as a new method facilitating the functional use of heat and 
opens a new possibility of spintronics [7] and spin caloritron- 
ics[8, 9]. 

Of particular interest is the SSE in the insulating mag- 
netic interface [4, 5]. The reason is that different from spin- 
dependent Seebeck effect, SSE allows heat to generate a pure 
flow of spin angular momentum, a flow of spins without elec- 
tron currents. This becomes obvious only after the observa- 
tion of the SSE through magnetic insulator-metal interfaces 
[4, 5]. Itinerant electrons are often problematic in the ther- 
mal design of devices, of which the issue can be avoided by 
the SSE in magnetic insulator interfaces without conducting 
electron currents. It allows us to construct efficient thermo- 
electric devices upon new principles [10] and to realize non- 
dissipative information and energy transfer in the absence of 
Joule heating [11, 12]. 

In this Letter, we predict new phenomena of the SSE across 
metal-insulating magnet interface: the rectification and nega- 
tive differential SSE. That is, reversing the thermal bias gives 
asymmetric spin currents and increasing thermal bias gives a 
decreasing spin current. The rectification of spin Peltier ef- 
fect (SPE) is also uncovered. We first demonstrate these inter- 
esting effects in a nanoscale magnetic interface system. Our 
simple model readily renders analytic interpretations and clear 
physical insights of their microscopic mechanism. We then 
turn to a macroscopic magnetic interface to demonstrate that 
the nontrivial rectification and negative differential SSEs are 
absent in such macroscopic system. Conditions to retain these 
intriguing properties are also discussed. The phenomena pre- 
dicted here are relevant for designing efficient spin/magnon 
diode and transistor, which could play crucial roles in spin- 
tronics [7], spin caloritronics [9] and magnonics [13], and 
could have potential applications in controlling energy and in- 
formation in low-dimensional nanodevices [14]. 

The system we are studying is schematically illustrated in 
Fig. 1, similar to the setup of longitudinal spin Seebeck exper- 
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FIG. 1 . Schematic illustration of the two-level system setup that de- 
scribes the nanoscale magnetic interface. 

iment [5]. We consider a nanoscale interface as in the situation 
of quantum point contacts. As such, the electron states near 
the interface can be simplified as impurity states or local states 
on a quantum dot: He = s^d^d^ + s^d^d^, where d^(d a ) is 
the creation (annihilation) operator of the local electron with 
spin a and energy e a at the interface. 

The left metallic lead is described in terms of the free 
electron gas: H L = T,k*( £ kcr - Ha)c\ a Ck*, O =t,l), 
with possibly different spin-dependent chemical potentials /j, a 
induced by spin accumulation [15, 16]. The spin voltage 
— fi^ can be measured by Hanle method [3] or be con- 
verted into an electric voltage through the inverse spin Hall 
effect [4, 5]. The electrons are freely exchanged between 
the left metal and the central quantum dot without spin flip: 
VL = T, k * t ka ci a d a + H.c. 

The right insulating magnetic lead is described by a Heisen- 
berg lattice: 

Hr = -J^\ S t S I + \ S tSj + SfSJ], (1) 

where 5^ is the raising (lowering) operator for the localized 
spin at site j, SJ is the spin operator of the z direction and 
J denotes the exchange coupling strength. It is convenient 
to map the spin operators into bosonic magnons by Holstein- 

Primakoff transformation [12, 17]: Sl~ = ^2 So — a^aj aj, 

S~ = a^2So — Ojdj, Sj = So — a^ctj, where So is the 
length of localized spins. Clearly, the creation (annihilation) 
of a local magnon a^(aj) at site j corresponds to the low- 
ering (rising) of the local spin component, i.e., excitation of 
magnons means that spins point less in z-direction and mag- 
netization goes down. At large spin limit or low temperatures 
((a^dj) <C 25q) we can approximate S~ w \/2So and 
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Sj~ ~ V^So a r Therefore, after a Fourier transform into 
the momentum space, the right insulating magnetic lead is ap- 
proximated by the free magnon gas: 

Hr « ^ ftw q a q a q + const., (2) 

q 

where the dispersion of uj q depends on the lattice details that 
are irrelevant to our present problem. 

Similar to Refs. [18-20] , we describe the electron-spin 
interaction at the nanoscale magnetic interface by the s-d ex- 
change coupling [21]: 

V R = -J2 Jj5|(4^-4^) + 5-4^ + 5+4dt], (3) 
q 

where J q denotes the effective exchange coupling at the in- 
terface, S~ ~ y/2So a J and « \^2Soa q . It is clear that 
the first term just splits the two local energy levels at interface 
so that its contribution can be absorbed into a renormalized 
e a (we keep the same notation) and we have the spin split- 
ting > generally. The second term S~dtdi describes 
the magnon emission process (into the right lead) associated 
with scattering of a spin-down electron to a spin-up electron 
at the interface. The last term S+d^d^ describes the reversed 
process that the spin-up electron in the dot absorbs a magnon 
from the right lead and excites to the upper level with spin-flip. 

Considering each magnon carries an angular momentum of 
—h (associated with a magnetic moment), the magnon current 
is equivalent to a spin-down current, which is then obtained by 

the Heisenberg equation of motion 1$ := \ ([ ? y R ] ^ 

or Is := %([Vr, ^2 q ^ q a q ]). From either definition, we can 
get the same spin (magnon) current: 

Notice that the expectation value (S~dldi) is just the com- 
plex conjugate of (S+dld^), as a next step, one has only to 
evaluate the equation of motion: 

j t (S+d\d t ) = 1 -{[H C + V R + H R , S+d\d t }) (5) 

= '-(e, - s t - fi w ,)(5+4d t > + ^,{[5+4^, 5"4dJ>. 

Note that the argument t of each operator is hidden 
for clarity. Following Ref. [22], in the steady state 
of large t limit, one can write (S+d^d^S q dLdi) = 

(S+S-)(4^)(V|> = 2 ^o(l + N R )f Li (l - f Lt ) where 
N R = [ e H/(^T R ) _ is me ma g non distribution 

function in the right lead, which acts as a thermal bath 
with temperature T R \ f La = [ e ^<r-^)/{k B T L ) + de _ 
notes the electron distribution function at the central dot, 
which is in equilibrium with the left lead at temperature 
T L . Similarly, (S-^S+djdt) = ( S q S q)( d i d l)( d \ d t) = 



2SoNRfL^(l — fLi)> Thus, by applying the Markov condi- 
tion in the steady state t oo, one obtains: (S+d]d^) = 
fcSoJg dTe^i-^-^WtfuO. - / Lt )(l + JV fl ) - 
/i/t(l — Il^Nr]. Substituting it into Eq. (4), we finally 
arrive at the expression of spin current: 
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is = ^r J [/^(i-/ Lt )(i+jv fl )-/ Lt (i-/^)jv fl ], (6) 

where Tj = 2tt ^ J^(^ - e t - fio; g ). 

This expression is familiar from the results of applying 
Fermi's golden rule to the definition Eq. (4). The first prod- 
uct — /z,t)(l + -^r) describes the relaxation rate of the 
occupied higher spin-down state flipping to the empty lower 
spin-up state with emitting a magnon with energy hu q into the 
right lead. The second product /l^ (1 — /li)Nr reversely de- 
scribes the excitation rate of the occupied lower spin-up state 
flipping to the empty higher spin-down state with absorbing a 
magnon from the right lead. The two spin-flip processes, ac- 
companied by the magnon emission/absorption, conserve not 
only the angular momentum, but also the energy, which is im- 
posed by the delta function 5(e± — — hw q ). In this way, 
only magnons with energy Huj q = e± — is able to trans- 
fer through the magnetic interface. Therefore, for the energy 
current carried by magnons, we can also obtain similarly: 
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/ =-^r J (£ 4 -£ t )[/ u (i-/ it )(i+Ar R )-/ it (i-/ L4 )Ar R ]. 

(7) 

Considering n u = /io± A/i s /2, T L , R = T ± AT/2, the 
thermal-spin transport coefficients are conventionally consid- 
ered in the linear response regime: expanding to the first or- 
der of spin voltage bias Afi s = fj,± — fif and thermal bias 
AT = T L - T R , which yields 

(ij = ( an G kT ° ) (at/To) ' (8) 

where Q = ^/ (s inh[^] -sinh[^] +sinh[^]) 
is the spin conductance, generated by the spin voltage differ- 
ence, k — Q{e] r — e^) 2 /T denotes the heat conductance, 
produced by the temperature bias. S = — A/i s /AT|/ s =o = 
(ei — s^)/Tq is the spin Seebeck coefficient, depicting the 
power of generating spin voltage by the temperature bias. 
II = Iq/Is\ at=o = £4, — £t * s me s P m P e lti er coefficient, 
depicting the power of heating or cooling carried by per unit 
spin current. It is seen clearly that we have II = STq so that 
the Onsager reciprocal relation QSTq = QU is fulfilled in the 
present system [23]. 

In the linear response regime, the spin Seebeck and Peltier 
coefficients are symmetric when reversing A/i s — > —Afi s and 
AT — » —AT. However, as we will see soon, when we go to 
the nonlinear response regime, the rectification of SPE and 
SSE will emerge. In some cases, we can even have the nega- 
tive differential SSE. 

We first examine the SPE which gives the magnonic heat 
current by the spin voltage. Figure 2 illustrates the rectifi- 
cation of SPE in our nanoscale interface system: the posi- 
tive spin voltage produces a large heat current, thus acting 
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FIG. 2. Rectification of spin Peltier effect, (a) Heat current Iq as a 
function of spin voltage A/i s with /j,^ = fio =b A/i s /2, for varying 
e t4 . no = 0, Tl = T R = T = 300 K. (b) Heat current as a 
function of spin voltage for varying To. e± — 50 meV, — 25 
meV, fio = 0. For all calculations in this work, we set Yj — 2.5 
meV, So = 16 [19]. 



as a good spin-heater or spin-cooler; but reversing the spin 
voltage gives a less heat current, thus acting as a bad spin- 
heater/cooler or even a spin-thermal insulator. Since Is is pro- 
portional to Iq, Is vs A /jl s will exhibit the similar behavior, 
thus acting as a spin current diode. 

The rectification of SPE can be understood in the follow- 
ing: when the spin voltage is positive (/i^ > /i^), the elec- 
tron prefers occupying the higher spin-down state and empty- 
ing the lower spin-up state. Thus, the transfer is dominated 
by the magnon emission accompanied by the relaxation of 
spin-down flipping to the spin-up state, so that we have Iq w 
^TjOi - e t )(l + N R ) from Eq. (7). When the spin volt- 
age is reversed, the electron prefers occupying the lower spin- 
up state and emptying the higher spin-down state. Then, the 
transfer is dominated by the magnon absorption accompanied 
by the excitation of spin-up flipping to the spin-down state, so 
that we have Iq « — ^-Tj(ei — s^)Nr with the negative 
sign denoting the reversed current direction. Therefore, the 
rectification ratio is obtained as (1+Nr) /Nr = exp ( £ ^~^ ) • 
This indicates that increasing the spin- splitting e± — or de- 
creasing the temperature T can enhance the rectification of 
SPE, which is consistent with Fig. 2 (a) and (b), respectively. 
Note that heat and spin currents prefer Afi s > (fi^ > /i^) is 
a consequence of the pre-defined up-magnetization direction 
of the insulating magnetic lead. Reversing the pre-defined 
magnetization direction (as well the central two levels) will 
reverse the rectification. 

We next examine the SSE which generates the spin current 
by the thermal bias. Figure 3 shows the rectification of SSE 
clearly that the thermal-induced spin currents are asymmetric 
with respect to reversing the temperature bias. In particular, 
when the two levels are either both above the chemical po- 
tential si > > no [see case (i) in Fig. 3(a)] or both be- 
low /io > £4. > [see case (iv) in Fig. 3(a)], we can even 
have negative differential SSE: when increasing the temper- 
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FIG. 3. Rectification and negative differential spin Seebeck effect, 
(a) Magnon spin current Is as a function of the normalized temper- 
ature bias AT/2T with Tl,r = T ± AT/2, for varying e t4 : (i) 
si 50 meV, e t 10 meV;' (ii) e± 50 meV, e t 0; (iii) e± = 50 
meV, £ t -25 meV; (iv) £4. -25 meV, e t -75 meV. Other 
parameters are fi^ = ^ = 0, To = 300 K. (b) Magnon spin cur- 
rent as a function of the normalized temperature bias, for varying To. 
si — 50 meV, — 25 meV, — = 0- 



ature difference Tr — Tl = —AT, instead of observing an 
increasing spin current as in linear response regime, we get 
a decreasing and even vanishing spin current from the right 
magnetic insulator to the left metallic lead. Figure 3(b) shows 
that the negative differential SSE also exists for a wild range 
of temperatures. 

The emergence of negative differential SSE can be reasoned 
as follows: when Tr > Tl, the thermal bias drives spin cur- 
rent from the right to the left. If near the linear response 
regime, it is natural to have a positive differential SSE that 
lowering Tl will increase the thermal bias which in turn in- 
creases the spin current. If we further decrease Tl with as- 
suming two levels are both above (below) the chemical po- 
tential, the two states will be both depleted (occupied), which 
in turn severely suppresses the magnon emission/absorption 
process that requires the concurrence of one occupied and one 
empty state. As a consequence, the conductance will decrease 
although the bias increases. When the conductance decreases 
faster than the increasing bias, the negative differential SSE 
emerges. 

Finally, we turn to examine the possibility of rectification 
and negative differential SPE and SSE in macroscopic mag- 
netic interfaces. In this situation, the two-level central system 
mimicking the nanoscale interfaces is discarded. As such, the 
left metal Hl interacts with the right magnetic insulator Hr 
directly, with the interfacial electron- spin interaction: 



H s d 



k,q 



(9) 



where the term —J q S*(c\^Ck^ — C/4C/4) is customarily ab- 
sorbed into Hl. Following similar procedures as above, 
from the definitions I s = ^[H sd ^ q a \ a q] and Iq = 
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\[Hsd, Eg tKUq%a q }, we obtain: 
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Io = 
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duF R {u:) / de PL (e)W(e,uj), (10) 

J — CO 

POO POO 

/ dcoF R (uj)tuv / dep L (e)W(e,w),(U) 

JO J-oo 



where W(e,o;) = f Li (e + ^)[1 - f Lt (e)][l + 
JVfl(M] - /Lt(e) [1 - hi(e + fiw)] jV fl (fio;) with = 
[ e (x- Mff )/(feBT L ) + !]-i m dN R (hu) = [e h ^/^TR) _ 

Note W has the same physical meaning as discussed after 
Eq. (6). The only difference is that we now have one inte- 
gral of magnon energy and one integral of electron energy for 
the macroscopic interface. F r (uj) contains the bulk magnon 
density of state (DOS) and the electron-spin couplings, which 
is reminiscent of Eliashberg function in the field of electron- 
phonon coupling [21]. Pl(z) denotes the bulk electron DOS 
in the left metal lead. 

It is customary and legitimate to take a constant bulk 
electron DOS Pl(z) ~ Pl(po) because for a good metal, 
the integral over de converges within a thermal energy 
JzbTl around the chemical potential [21]. Then, by apply- 
ing the equality j defLi(s + hu)[l — fLt( £ )] — — 

Ap s ), we finally arrive at the formulas 

Is= 2SoPl(po) J~ fa FR{ „ ){hu , _ A/Ma) 

x[N L (ftLj-Ap s )-N R (ftu)]- (12) 
Iq = 2SoP ^ o) J°° duF R (u)hw(hu - Ap s ) 

x[N L (hu-Ap s )-N R (hw)}. (13) 

Clearly, in the absence of thermal bias Tl = T R , the rectifica- 
tion of SPE exists, i.e., the heat current Iq is asymmetric un- 
der reversing the spin voltage Ap s — » —Ap s . Similarly, Is is 
also asymmetric under revering the spin voltage. Moreover, it 
is straightforward to verify that Oa^^s and <9a/i s ^q are both 
positive. Therefore the negative differential SPE is absent. 
These observations are the same as above in the nanoscale 
magnetic interface. 

For SSE, we take Ap s = and Tl ^ T R . In this way, Eqs. 
(12, 13) reduce to Landauer-type formulas 

j^ ^OpMp ^^ [Nl{uj) _ Nr{uj)] . (14) 

_2S p L (p Q ) 



h 



7 dujF R (uj)(hw) 2 [N L (oj)-N R (uj)}.(15) 
Jo 



Since the temperature dependence only manifests in Bose dis- 
tributions Nl and N R , one can readily verify that for this case 
of bulk macroscopic magnetic interfaces we can never have 
the rectification and negative differential SSE. 

One may wonder how is the conclusion of the macroscopic 
magnetic interface so different from that of the nanoscale in- 
terface. From the above analysis, one can find that the con- 
stant bulk electron DOS is the crucial assumption that leads 



to the Landauer-type Eqs. (14, 15) from Eqs. (10, 11). In 
low-dimensional nanoscale systems, the DOS should not be a 
smooth function but vary strongly on energy. Therefore, when 
the magnetic interface scales down to the nanocale, we have 
to keep the electron DOS Pl(z) inside the energy integral. In 
this way, we will not have the Landauer-type equations but in- 
stead retain the intriguing properties of rectification and neg- 
ative differential SSE. 

For example, considering a common Lorentzian-type DOS 
where the half-width T is small, the 



PL(G) — k ( £ _ £o )2 +F 2 

resonant peak of DOS becomes sharp at sq, so that Eqs. (10, 
11) reduce to 



Is = 
Io = 



2So 

h 
h 



(16) 



POO 

/ dwF R (u;)W(eo,u>), 
Jo 

nOO 

/ duF R (u)hujW(eo,uj). (17) 
Jo 



These results are familiar from Eqs. (6, 7). The mere differ- 
ence is that the currents now have an integral over the magnon 
spectrum, instead of the resonant transfer of a single-mode 
magnon. It is straightforward to verify that Eqs. (16, 17) have 
similar behaviors of the rectification and negative differential 
SSE as Eqs. (6, 7). 

Generally, the Eliashberg-like function F R (u) also has 
electronic energy dependence [21]. By considering non- 
smooth, strong energy dependent F R (e,uj) in some designed 
magnetic interfacial systems, we expect to have the similar 
rectification and negative differential transport effects. 

In summary, we have studied the nonequilibrium spin- 
thermal transport across a metal-magnetic insulator interface. 
The conjugate-converted thermal-spin transport is assisted by 
the exchange interaction at the interface, between conduction 
electrons in the metal lead and localized spins in the insulating 
magnet lead. We have predicted the rectification and negative 
differential SSE as well as the rectification of SPE, and re- 
solved their microscopic mechanism. Since the heat current 
is proportional to the magnon spin current, we also possess 
the rectification of spin current driven by spin voltage, and 
the rectification and negative differential thermal conductance 
driven by the temperature bias, as in dielectric phononics [14]. 

It is also readily to generalize our findings to the interfaces 
of metal-magnetic metal/semiconductor. Since recent studies 
imply the important role of phonon-drag in SSE [6, 24-27], 
taking account of the effect of nonequilibrium phonons on rec- 
tification and negative differential SSE would be an interest- 
ing topic. By integrating the phononics [14] with spintronics 
[7], spin caloritronics [9] and magnonics [13], we expect to 
possess more opportunities to achieve smart control of energy 
and information in low-dimensional nanodevices. 
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